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Abstract 

In this paper we consider a simple Markov chain for bipartite graphs 
with given degree sequence on n vertices. We show that the mixing time 
of this Markov chain is bounded above by a polynomial in n in case of 
semi-regular degree sequence. The novelty of our approach lays in the 
construction of the canonical paths in Sinclair's method. 

1 Introduction 

The degree sequence, d{G), of a graph G is the non- increasing sequence of its 
vertex degrees. A sequence d = (c?i, . . . , d„) is graphical iff d{G) = d for some 
simple graph G, and G is a graphical realization of d. 

Already at the beginning of the systematic graph theoretical research (late 
fifties and early sixties) there were serious efforts to decide whether a non- 
increasing sequence is graphical. Erdos and Gallai (1960, [4]) gave a necessary 
and sufficient condition, while Havel (1955, [6 ) and Hakimi (1962, J]) inde- 
pendently developed a greedy algorithm to built a graphical realization if there 
exists any. (For more details see for example [TU].) 

Generating some (or all possible) graphs realizing a given degree sequence or 
finding a typical one among the different realizations is an ubiquitous problem 
in network modeling, ranging from social sciences to chemical compounds and 
biochemical reaction networks in the cell. (See for example the book [13] for a 
detailed analysis, or the paper |10j for a short explanation.) 

When the number of different realizations is small, then the uniform sam- 
pling of the different realizations can be carried out by generating all possible 
ones and choosing among them uniformly. 

However in cases where there are many different realizations this approach 
can not work. In these cases some stochastic processes can provide solutions. 
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Here we mention only one of the preceding results: Molloy and Reed (1995, [T^ ) 
applied the configuration model (BoUobas (1980, [IJ) for the problem. They 
successfully used it to generate random graphs with given degree sequences 
where the degrees are (universally) bounded. Although it is well known, that 
this method is computationally infeasible in case of general, unbounded degree 
sequences. 

A completely different method was proposed by Kannan, Tetali and Vempala 
(1999, [5]), which is based on the powerful Metropolis-Hastings algorithm: some 
local transformations (these are the well-known swaps) generate a random walk 
on the family of all realizations. They conjectured that this process is fast mixing 
i.e. starting from an arbitrary realization of the degree sequence the process 
reaches in a completely random realization in reasonable (i.e. polynomial) time. 
They succeeded to prove it for bipartite regular graphs. Their conjecture was 
proved for arbitrary regular graphs by Cooper, Dyer and Greenhill (2007, [5]). 

The goal of this paper was to attack Kannan, Tetali and Vempala's conjec- 
ture for arbitrary bipartite degree sequences. We obtained the following result, 
which is one further step toward the general Kannan, Tetali and Vempala's con- 
jecture. A bipartite graph is called semi-regular if in one class all the vertices 
have the same degree, while the degrees in the other class are arbitrary. 

Theorem 1.1. The Markov process - defined by Kannan, Tetali and Vempala 
- is fast mixing on bipartite semi-regular degree sequences. 

This paper is rather long therefore at first we give a short sketch what is going 
to happen later. 

2 Outline of the paper 

Assume that G and G" are two realizations of the same bipartite degree sequence 
(see formula (|3.ip ). The realizations are given by the mean of adjacency matrices 
Mg of the graph G. As Ryser proved ([Ilj, see Theorem 13. 2p they can be 
transformed to each other with a sequence of swaps (for exact description see 
Definition 13.11 ). In this paper we give a short, transparent proof of this fact 
(using the language of graphs instead of the language of 0-1 matrices). 

Since the inverse transformation of a swap is also a swap, therefore the 
Markov graph G (see Section 0]) - which consists of all possible realizations of 
a fixed bipartite degree sequence - is connected. One can easily prove that G's 
diameter is at most the sum of the elements in the degree sequence (see the 
Proof of Theorem l3.2[) . Equipping with the appropriate transition probabilities 
(see formula (|4.H) ). this depicts a Markov chain on the set of all realizations. It 
is easy to see that this Markov process is reversible with the uniform distribu- 
tion as the globally stable stationary distribution. Therefore an application of 
the Metropolis - Hastings algorithm will converge to some completely random 
element, independently from its starting point. The only open problem here is 
its mixing time: how fast a random sample can be reached. 
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In their papers, Kannan, Tetali and Vempala (1999, [5]) (and Cooper, Dyer 
and Greenhill, 2007, [2]) used this method to sample almost uniformly all pos- 
sible realizations of bipartite regular (general regular) degree sequences, resp. 
Both approaches are based on the seminal Sinclair's method (Section [5]). 

For the successful application of this method you have to construct a set of 
canonical paths, which satisfies a long list of properties (see Section [S]). The key 
property is the following: if X and Y are two realizations of the degree sequence, 
then for each canonical path P (leading from X to Y) and for all vertices Z in 
P there exists a realization U such that these four realizations together satisfy 
the condition: 



where C is a fixed constant (here d{M, M') denotes the number of different 
entries in two integer matrices). The successful construction of the canonical 
path system provides the necessary input of the Sinclair's method, which in 
turns proves the fast mixing behavior of the Markov process. 

In the canonical path construction process, for all pairs X and Y of possible 
realizations we have to define a collection of swap sequences leading from X to 
Y (and, automatically, from Y to X). For that end one can take the symmetric 
difference of the edges of the two realizations, then it can be uniquely decom- 
posed into cycles (with edges from X and Y altcrnatingly, see the list " Outline of 
the construction of the path system" in SectionlHl). In this point both the Tetali 
and Vempala and the Cooper, Dyer and Greenhill methods apply a uniform 
"folding" process on those cycles to determine the swap sequences. 

The main contribution of our paper is the introduction of the friendly path 
method to construct suitable canonical path systems. Our method takes into 
consideration the environment of the alternating cycles. For semi-regular bi- 
partite degree sequences one can prove that our construction satisfies condition 
(j2.1l) therefore the Markov process under review is fast mixing. However when 
the edge sets of the bipartite graphs X and Y differ precisely in one (alternat- 
ing) cycle, then the friendly path method always provides a suitable canonical 
path system. (See Theorem 18. ip . 

To finish this Section we also want to mention that any graph can be repre- 
sented as a semi-regular bipartite graph: here one class represents the vertices, 
while the other class describes the (labeled) edges. Using this representation for 
sampling randomly the original graphs faces the problem that different graphs 
will appear with different frequencies. An appropriate weighted Metropolis- 
Hastings algorithm may provide a solution. 

3 Basic definitions and preliminaries 

Let G — ([/, V; E) be a simple bipartite graph (no parallel edges) with vertex 
classes U — {ui, . . . ,u{\,V — {wi, . . . , w^}. The (bipartite) degree sequence of 
G, bg{G) is defined as follows: 



d{Mx + My - Mz,Mu) <C 



(2.1) 



bd{G) = [{diu,),...,diuk)),{divi),...,d{vi))), 



(3.1) 
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where the vertices are ordered such that both sequences are non-increasing. 
From now on when we say "degree sequence" of a bipartite graph, we will 
always mean the bipartite degree sequence. 

A pair (a, b) of sequences is a (bipartite) graphical sequence (BGS for short) 
if (a, b) = bd (G) for some simple bipartite graph G, while the graph G is a 
(graphical) realization of (a, b). 

Next we define the swaps, our basic operation on bipartite graphs. 

Definition 3.1. Let G — [U, V, E) be a bipartite graph, ui, U2 G U, vi,V2, G V , 
such that induced subgraph G[mi, U2, ; wi, W2] is a 1-factor, (i.e. (ui,Vj), (u2, ws-j) G 
E, but {ux,vz-j),{u2,Vj) ^ E for some j.) Then we say that the swap 
on (ui, U2; fi, wi) is allowed, and it transforms the graph G into a graph 
G' = {U,V,E') by replacing the edges t;^), (u2, ua-j) by edges (ui,U3_j) 
and (w2, Vj), i.e. 

E' = E\{{ui,Vj), (u2, W3-i)} U {(wi, W3-j), (u2,Vj)}. (3.2) 

So a swap transforms one realization of the BGS to another (bipartite graph) 
realization of the same BGS. The following proposition is a classical result of 
Ryser (1957, [H]). 

Theorem 3.2 (Ryser). Let Gi = [U, V] Ei) and G2 = ([/, V] E2) be two realiza- 
tions of the same BGS. Then there exists a sequence 0/ swaps which transform 
Gi into G2 through different realizations of the same BGS. 

Ryser's result used the language of - 1 matrices. Here, to make the paper self 
contained, we give a short proof, using the notion of swaps. The proof is based 
on a well known observation of Havel and Hakimi ( [H [5] ) : 

Lemma 3.3 (Havel and Hakimi). Let G = {U,V; E) be a simple bipartite graph, 
and assume that d(u') < d{u), furthermore u'v £ E and uv ^ E. Then there 
exists a vertex v' such that the swap on (u, u'; v, v') is allowed, and so it produces 
a bipartite graph G' from G such that Tgi{v) — {Tg{v) \ {u'}) U {u}. 

Proof: By the pigeonhole principle there exists a vertex v' ^ v such that uv' £ E 
and ulv' ^ E. So swap defined on vertices (u, w'; v, v') is allowed. □ 

We say that the previous operation is pushing up the neighbors of vertex v. 
Applying the pushing up operation d times we obtain the following push up 
lemma. 

Lemma 3.4 (Havel and Hakimi). If G ~ {U,V;E) is a simple bipartite graph, 
d{ui) > d{u2) > • • • > d{uk) and v € V , d = d{v). Then there is a sequence 
S of d many swaps which transforms G into a graph G' such that Fg' (w) = 

{iti, . . .,Ud}. 

This pushing-up lemma also suggests (and proves the correctness of) a greedy 
algorithm to construct a concrete realization of a BGS (a, b) . 

Proof of Theorem 13.21 We prove the following stronger statement: 
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there exists a sequence of 2e swaps which transform Gi into G2, where e 
is the number of edges of Gi 

by induction on e. Assume that {^b) holds for e' < e. We can assume that 
d{ui) > d{u2) > ■ ■ ■ > d{uk))- By the Havel- Hakimi Push-up Lemma there is a 
sequence Ti (T2) oi d = d{vi) many swaps which transforms Gi (G2) into a G'l 
{G'2) such that Ta^ivi) = {ui, ...,Ud} iTc'^{vi) = {ui, . . .,Ud}). 

We consider the bipartite graphs G'{ = G[ \ {vi} and G2 = G2 \ {wi}, i-e. we 
remove the vertex vi and all the edges connected to vi. Since hd{G'l) — ^^(Gj) 
and the number of edges of G'l is e — d, by the inductive assumption there is a 
sequence T of 2(e — 0?) many swaps which transforms G'/ into Gj- 

Now observe that if a swap transforms H into iJ', then the "inverse swap" 
(choosing the same four vertices, and changing back the edges) transforms H' 
into H . So the swap sequence T2 has an inverse which transform G2 into G2. 
Hence the sequence T1TT2 is the required swap sequence: it transforms Gi into 
G2 and its length is at most d -f 2(e — d) -\- d = 2e. □ 



pala's lead - we define a Markov chain (G, V) in the following way. G is a graph, 
the vertex set V{G) of the graph G consists of all possible realizations of our 
BGS, while the edges represent the possible swap operations: two realizations 
are connected if there is a swap operation which transforms one realization into 
the other (and, recall, the inverse swap transforms the second one to the first 
one as well). 

We use the convention that upper case letters X, Y and Z stands for vertices 

of y(G). 

The graph G clearly may have exponentially many vertices (that many dif- 
ferent realizations of the degree sequence). However, by the statement (^) (in 
the proof of Theorem 13. 2|) . its diameter is always relatively small: 

Corollary 4.1. The swap distance of any two realizations is at most 2e, where 
e is the number of edges. 

V denotes the set of the transition probabilities, which are defined as follows: 
if the current realization (status of the process) is G then we choose uniformly 
two- two vertices iti, 1*2; vi, f2 from classes U and V respectively and perform 
the swap if it is possible and move to G'. Otherwise we do not perform a 
move. The swap moving from G to G' is unique, therefore the probability of 
this transformation (the jumping probability from G to G' ^ G) is: 



The probability of transforming G' to G in any time is just the same. The 
transition probabilities are time and edge independent and they are also sym- 
metric. Therefore V can be imagined as a symmetric matrix, where all off 




Prob(G^>G') :=r(G'|G) 



1 



(4.1) 
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main-diagonal, non-zero elements are the same, while the entries in the main 
diagonal are non-zero, but (probably) different values. 

As we observed, the graph G is connected, therefore the Markov process 
is irreducible. Furthermore since in each realization (state) our random walk 
can stay in that state with positive probability (except for the unique degree 
sequence ((1, 1), (1, 1))), therefore our Markov process is clearly aperiodic. 

Finally since, as we saw, the jumping probabilities are symmetric T{G\G') = 
T{G'\G), therefore our Markov process is reversible with the uniform distribu- 
tion as the globally stable stationary distribution. 

Since our Markov process is determined by the graph G, sometimes G will 
be mentioned as the Markov graph. 

5 The Sinclair's Method 

To start with we recall some definitions and notations from the literature. Let 
denote P* the tth power of the transition probability matrix and define 

A,(i) :=imax|P*(y,x)-l/iV|, 

and 

T^(e) := min{A2;(t') < s for all t' > t}. 
Our Markov chain is said to be fast mixing iff 

T.(e) <o(poly(/og(iV/e))). 
Consider the different eigenvalues of P in decreasing order: 

1 = Ai > A2 ■ • • > Am- 

The relaxation time Trei is defined as 

_ 1 

- T^- 

The following result was proved by Diaconis and Strook in 1991: 

Theorem 5.1 (Diaconis and Strook, |3]). Tx{e) < Trei ■ poly(ln£^^ + In A^). □ 

So one way to prove rapid convergence is to find a polynomial upper bound of 
Trei and we need fast convergence of the method to the stationary distribution 
otherwise the method cannot be used in practice. 

Kannan, Tetali and Vempala in ^ could prove that the relaxation time of 
the Markov chain {G,V) is a polynomial function of the size n := 2k of (a, b) 
if it is a regular bipartite degree sequence. Here we extend their proof to show 
the fast convergence of the process for the semi-regular bipartite case. 

There are several different methods to prove fast convergence, here we use - 
similarly to - Sinclair's canonical path method ([15j). 
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Theorem 5.2. Let W be a graph whose vertices represent the possible states of 
a time reversible finite state Markov chain A4, and where uv G ^^(H) iff the 
jumping probabilities in A4 satisfy T{u\v)T{v\u) ^ 0. For all x ^ y E V(M.) 
denote T^.y a set of paths in M connecting x and y together with a probability 
distribution n^^y on T^^y. Furthermore let 

r := U r.,, 

x^yeV{U} 

where the elements ofT are called canonical paths. We also assume that there 
is a stationary distribution tt on the vertices V(H). Finally let 

Kr '■= max 7r(a;)7r(w)7rj^ ,,(7) — — — - — ; — -. (5.1) 

x,yiiV{U) (u.,z)G£;(H): V I y V ; 

76r^_j : ee7 (■u;,z)e7 

Then 

r,ei(X) < (5.2) 
holds. □ 



We are going to apply Theorem 15.21 for our Markov chain (G,!^). Using the 
notation |V(G)| :— N, the (uniform) stationary distribution has the value 
7r(X) = 1/N for each vertex X G V{G). Furthermore each jumping proba- 
bility has the property T{x\y) > 1/n'^. So if we can design a canonical path 
system such that each path is shorter then an appropriate poly(n) function, 
then simplifying inequality ()5.ip we can turn inequality ()5.2p to the form: 



Trol < 



/ \ 

poly(n) , . 

max } T^xAl) 

x.,y&{H) . 



N 



(5.3) 



So our task is to find a canonical path system {Tx.y : x,y £ G} such that 

• each path is shorter then an appropriate poly(n) function, 

• in (|5.3p we have a polynomial upper bound for the right side. 



6 Canonical paths - general considerations 

Our construction method for canonical paths commences on the trail of Kannan, 
Tetah and Vempala ([5]), and Cooper, Dyer and Greenhill ([2])- 

li X,Y e V{G) let E{X A Y) be the symmetric difference of E{X) and 
E{Y), set E{X-Y) = E{X)\E(Y), and E{Y ~ X) ^ E{Y)\E{X). The edges 
in E{X — Y) are the X-edges and the edges in E{Y — X) are the Y -edges. Due 
to the fact, that the two degree sequences are equal, for every vertex u £U yjV 
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the number of X-edges adjacent to u is the same as the number of F-edges 
adjacent to u. 

Before we describe the construction of our canonical path system we fix some 
definitions: 

Definition 6.1. In a simple graph a circuit is a sequence of pairwise disjoint 
edges ei, . . . , e^, where each pair of edges ej, Cj+i (the summation goes by modulo 
£) share a common end point, furthermore the common end points of ei-i,ei 
and Ci, Bi+i are different. If these common end points are pairwise distinct, then 
this circuit is a cycle. 

We start our construction process by defining 

S{u; {E{X-Y),E{Y-X))) = 

= {/ U /-^ I / : E{X - Y){u) E{Y - X){u) is a bijection} 

for each vertex u. The set <S(m; {E{X — Y),E{Y — X))) is the set of all possible 
bijections from the X-edges adjacent to u into the F-edges adjacent to u. We 
also define 

§(X, Y) =llS{u; {E{X - Y),E{Y - X))) . 

u 

and 

§ = |J{§(X, Y):X,Y& V{G.)}. 

For s G § let 

A, = |J{dom(s(«;)) : w G f/ U F}. 
If s G S(X, Y) then clearly 

l^s=E{X/\Y). (6.1) 

If M and M' are mxm' matrices then let 0(M, M') be the number of non-zero 
elements in M — M'. For G G G let Mq be the bipartite adjacency matrix of 
G. For X,Y,Z gG write Mx,y,z = Mx + My - Mz- 

Outline of the construction of the path system. 

(A) To each s G S assign a circuit decomposition 

of Ag as follows: 

(a) Consider the graph Q ~ ( Ag , £) , where 

((u, v), {u'v')) G £ iS u = u' and s{u){u, v) = {u, v'). 

(b) 5 is a 2-regular graph because (u, v) G A^ is adjacent to s{u){u, v) and 
s{v){u, v), so it is the union of vertex disjoint cycles {Wf, . . . , }. 
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(c) The cycles of Q are circuits in E{X AY), so we set {VTf , . . . , W^^} as 
our desired decomposition of As into circuits. 

(B) For each X,Y eG and s e S{X, Y) 

A, ^E{X AY), 

by (|6.ip and the circuits 

form a decomposition of E{X A Y) into alternating circuits. 

(C) Decompose each Wf into alternating cycles 

1 j'~'2 T ■ ■ ^ ^l^^i' 

(D) Let 

J , 02 , . . . , . 

be the short hand notation of the cycle decomposition 

^1 1^2 I • ■ • I '-'^3 1 1 ^1 ' *-^2 7 • ■ • 7 '-^^g 2 ' ■ ■ ■ ' 1 ' 2 ' ■ ■ ■ ' fca 

of E{X AY). 

(E) For each cycle C in this decomposition we fix a walk-around of the cycle. 
If a,b £ C let [a, 6]c be the walk from a to 6 in C according to the fixed 
walk-around. 

(F) Define the path T{X, Y, s) 

X — Go, Gi, . . . , Gni 7 Gni+i, . ■ . , G„2, . . . , G„^^ — Y 
in G from X to Y such that 

(a) < C ■ n"^. 

(b) £;(G„,^J = {E{Gn,) u n £;(g,)) \ {e{x) n £;(g,))7 

(c) if rii < j < TT-i+i then there are two vertices a and b in G„. such that 

\E{Gj)AF\ <ni, 

where 

F = (i?(G„J u ([a,6]c, nii;(r)) \ (K&]c, ni?(x))), 

(d) for each i there is Hi <E G such that 
where fi,- are fixed "small" natural numbers. 
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(G) For 7 e Txx let 

, ^ \{sGS{X,Y):TiX,Y,s)=^}\ 
^.Al) = ^^^^ (6.2) 

Key Lemma 6.2. If we can assign paths 

{T{x, r, s) : s e r), r e g} 

according to (A)-(F) t/ien for each e G _E(G) 

X^I^^.^-It) : e e 7 e rjf,y : a;, y e G} < poly(n) • N, (6.3) 
and so our Markov chain is fast mixing. 
To prove this statement we need some preparations: 

Lemma 6.3. There is function ^' and a parameter set B such that B has poly(n) 
elements, and for each X,Y <E V{G), for each 7 = T{X,Y,s) G T{X,Y) and 
for each edge e ~ {Z, Z') from 7 there is B eM such that 

* (Z, Mx^Y,z, s, b) = {X, Y, 7). (6.4) 

Proof: By definition s determines the circuit decomposition {Wf, . . . , ) of 
As, and this circuit decomposition in turn determines a canonical alternating 
cycle decomposition (Cf , . . . , C^J by (A)-(D). 

Assume that Z = Gj in the path T{X, Y, s) such that Ui < j < n^+i. Put 



\J{E{a,)nE{z)} u[[a,b]c,nE{Z) 

i' <i 

u[[b,a]cAE{Z)) U [j{Eia,)\E{Z)} 



i'>i 

Then, by (F)(b)-(c) 

\{E{X)\E{Y)) ^F\ < n. 

So if the parameter B contains j, a, b and the set {E{X) \ E{Y)) A F then we 
can determine E{X)\E{Y). Clearly the size of the parameter set is polynomial. 
Since Z and Mx,y,z determine E{X) r\E{Y), we can compute E{X). Similarly 
we can compute E{Y). Then path 7 = T(X, F, s) is uniquely determined by X, 
Y and s. □ 

Proof of the Key Lemma: Using these observations, we can rearrange HS 
of (16.31) as follows: 



^{7r;,.^(7) : X, r G V{G) and 7 G Tx,y : e G 7} = 
{^^•^i^) 3^^; 3s; 3B3Z ee 

s.t. vI/(z;A#;s,i?) -(X,r,7)} (6.5) 



10 



(It is worth to mention that it may happen, that - as we already mentioned it - 
several s give rise to the same path 7, but ttx.y (7) is still taken into consideration 
only once in the sum of LHS.) Let 

m = ^Mx.Y^z ■■ Z e T{X, Y, s) for some X, F e G and se S{X, 

By (F)(d) each M is in a small radius neighborhood of some ?7 e G, so 

\m\<poly{n)-\ViG)\=po\yin)-N. 
For M e an and B G B write 

It is enough to prove that 

r(M,S) <poly(n). 

So fix M and B. Let 

X(M,S) = {(X,y,s) : *(z,M,s,b) =(X,r,7)Aee7}. 

We can assume X 7^ 0. As we observed, 

E{X AY) = As. 
E(X A Y) has a degree sequence {2di, . . . , 2dh). Put 

h 

1 

Since = ioi each {X,Y,s) £ X, expression (|6.2p can be written as 

|{g£S(X,y)) : TiX,Y;.s)=^}\ 

T^x.Yh) = • 

Then 



T(M, B) 



However, the RHS of the previous formula is just |X|/iA- So it is enough to 
prove that 

|X(Af,B)| <poly(n)-U. 
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Since e, M, s and B determine X and Y it is enough to prove that 

|s : e G (X,r,s) £ X(M,B)| < poly(?i) -iA- (6.6) 

Assume that {X, Y, s) G X{M, B). Assume that Z = Gj and Ui < j < rij+i 

Then we can observe that s is "almost" in §(A n E{Z),A \ E{Z)). More 
precisely there is s' G §(A n A \ ^^(Z)) such that 

(Us) A (Us') I < 6* (6.7) 

for some constant 9. Indeed, by (F)(c), apart from a few edges, for i' ^ i the 
cycles Ci' alternates between Ani?(Z) and A\i?(Z)) by (F)(c). Moreover both 
[a, 5]ci and [&, ajc; are "almost alternate. This proves (|6.7I) which implies (j6.6p 
because |S(An£;(Z), A \£;(Z))| =iA. □ 

We try to carry out the plan we just described. So; 

• Fix X, r e G. 

• Pick s e §(A:,y). 

• s gives an alternating cycle decomposition 

Co,Ci,...,C,. (6.8) 

oiE{X l\Y). 
We want to define a path 

X — Go , . . . , Gi , . . . , G„i = Y 
from X into F in G - denoted by T{X, Y, s) - such that 

(i) the length of this path is < C ■ n? , 

(ii) for some increasing indices < ni < 712 < . . . we have G„. — Hi, where 



E{H.,)=E{X)a([_]e{C,)\ 

\i'<i / 



So we have certain "fixed points" of our path T{X,Y,s), and this observation 
reduces our task to the following: 

• for each i < I construct the path 

TT /-^/ r^i rr 

between G„; and G„.^j such that m' is < \Ci\ and (F)(d) holds, i.e. for 
each j there is Kj G G such that 0(Afx,y.G' 1 Kj) ^ ^3- 
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From now on we work in that construction. To simplify the notation we write 
G = Hi and G" = We know that the symmetric difference of G and G' is 

just the cycle Now we are in the following situation: 

Generic situation - construction a path along a cycle 

(i) X, Y, G, G' e G. 

(ii) The symmetric difference of E{G) and E{G') is a cycle C. 

(iii) the symmetric differences E{XAG), E{GAG') and E{G' AY) are pairwise 
disjoint. 

Construct a path 

G^Go,...,Gm^G' (6.9) 
in the Markov graph G such that 

(a) TO < c - |C|, 

(b) for each j there is Kj G G such that 'i){Mx,Y.Gj , ^k,) < ^^4- 

We will carry out this construction in the next chapters. The burden of such a 
construction is to meet requirement (b). In [9J and in [2 the regularity of the 
realizations was used. 

The "friendly path method". 

In the next chapters we describe a new general method based on the notion of 
friendly paths (see Definition 17. 3p to construct the paths T{X,Y, s). 

The novelty of our friendly path method can be summarized as follows: 

• if our bipartite degree sequence is semi-regular then the paths T(X, Y, s) 
satisfy (b) 

• if our bipartite degree sequence is arbitrary, then T{X, Y, s) satisfies (b) 
provided the symmetric difference of X and y is a cycle. (You can not 
get such paths using the methods of [S] and [2].) 

Originally we conjectured, that our friendly path method always produces paths 
which satisfy (b). However we were unable to prove it, and now we think that 
essentially new ideas are needed to prove the case of general bipartite degree 
sequences. 

7 Canonical paths - along a cycle 

In this section we describe the construction a canonical path along an alternating 
cycle C. 

In what follows we will imagine our cycles as convex polygons in the plane, 
and we will denote the vertices of any particular cycle of 2£ edges with ui,vi,U2, 
V2, ■ ■ . ,Ui,ve The edges of the cycle are (wi, tji), (ui, U2), . . . , {ui, Vi), (ve, ui) and 



13 



they belong alternately to X and Y. All the other (possible, but not necessarily 
existing) edges among vertices of a particular cycle are the chords. (With other 
words we will use the notion of chord if we want to emphases that we do not 
know whether the two vertices form an edge or not in the current graph.) A 
chord is a shortest one, if in one direction there are only two vertices (that is 
three edges) of the cycle between its end points. The middle edge of this three 
is the root of the chord. 

When we start the process we have realization G, when we finish it we have 
realization G' . These two have almost the same edge set, except that the X- 
edges in C belong to G while G' contains the F-edgcs. W.l.o.g. we may assume 
that the vertices are ui,vi,U2,V2, ■ . ■ ,Ui,Vi and {ui,vi) is an edge in G while 
{vi,U2) belongs to G'. We are going to construct now a sequence of graphical 
realizations between G and G' such that any two consecutive elements in this 
sequence differ from each other in one swap operation. The general element of 
this sequence will be denoted by Z. 

We have to control which graphs belong to this sequence. For that purpose 
we assigned a matrix M to each graph Z. If G is a vertex in G then denoted Mq 
the adjacency matrix of the bipartite realization G where the columns and rows 
are indexed by the vertices of U and V resp., but the columns are numbered 
from left to right while the lines are numbered from bottom to the top. Then 
let 

M{X + Y - Z) = Mx + My - Mz. 

By definition each entry of an adjacency matrix is or 1. Therefore only 
— 1,0, 1,2 can be the entries of M. An entry is —1 if the edge is missing from 
both X and Y but it exists in Z. This is 2 if the edge is missing from Z but 
exists in both X and Y. It is 1 if the edge exists in all three graphs {X, Y, Z) 
or it is there only in one of X and Y but not in Z. Finally it is if the edge is 
missing from all three graphs, or the edge exists in exactly one of X and Y and 
in Z. (Therefore if a chord denotes an existing edge in exactly one of X and Y 
then the entry corresponding to this chord is always or 1.) 

Observation 7.1. The raw and column sums of M{X + Y — Z) are the same 
as raw and column sums in Mx (or My or Mz). □ 

We need some more notions: 

Definition 7.2. The type of a chord is 1 or depending whether it is an edge 
in G or - what is the same - in G'. A chord / = {ua, vp) is a cousin of a chord 
e = {ue, vs), if < 6 — a < 1 and 0<e — /3<1. A chord e is friendly if at 
least one of its cousins has the same type as e itself. 

Since the notion of the cousins is not a very natural one we want to illustrate it. 

Therefore now we introduce our main tool what we will use later in this paper to 
illustrate different procedures in our current realizations. This tool F(G+G' +Z) 
(or Fz for short, and we use F for "figure") is a matrix again, similar to the 
matrix Mq + Mc + Mz. This always illustrate the subgraph in Z spanned by 
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the vertices of one of the alternating (therefore of an even cardinahty) cycle 
Ci. The positions (1, 1), . . . , {£,£) form the main-diagonal while the positions 
right above the main-diagonal as well as the rightmost bottom one (these are 
(1, 2), (2, 3), ...,(€- 1, £) finally {£, 1)) form the small-diagonal. 

The matrix Fz^s entries in the main-diagonal and the entries in the small- 
diagonal come from the adjacency matrix of Mz- The other entries come from 
the corresponding entries of Mq+Mg' +Mz- In that way in the main- and small- 
diagonal's elements are or 1 while the others (the off-diagonal entries) can 
be 0, 1, 2, 3. (There is an easy algorithm to construct Fz from the corresponding 
M{G+G' — Z) and vice versa (please recognize that here we use G and G" instead 
of X and Y): In the main-diagonal and in the small-diagonal the zeros and ones 
must be interchanged to the opposite ones. Outside of these diagonal entries 
-1, 0, 1, 2 of M{G + G' - Z) become 1, 0, 3, 2 in Fz-) Since in this ahernating 
cycle G and G" contain the same chords therefore the off-diagonal elements in 
Fz are odd when the edge exists in the actual Z and even otherwise. When 
Z = G then the main-diagonal entries are 1 while the small-diagonal elements 
are 0. This matrix Fz will be used in our illustrating figures. 

Now Figure [1] illustrates the cousins of the chord {uq,V2) in the initial Z. 
(They are (ui, wg), {ui,vr), {u2, vq), finally (m2, w?) and let's recall that the word 
chord indicates that the definition does not depend on the actual existence or 
non-existence of that edge.) 



Figure 1: An edge and its cousins 
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Definition 7.3. A sequence of pairwise distinct chords ei, . . . , is a friendly 
path if 

(i) each chord is friendly, 

(ii) each Ch has preciously one common endpoint with eh+i while the other 
end points are in distance 2, 

(iii) finally ei and Cj are shortest chords and their roots have different types. 
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A friendly path can be quite complicated, and it is important to remark that 
such a friendly path is NOT a path in a particular graph. The name is justified 
by the image of the friendly path in the illustration of Fz what we are going to 
show next. (It shows the path itself, but it does not show why the individual 
elements of the path are friendly - see Figure [51) This figure is not for illustration 
only: whenever we consider a friendly path we always work on the matrix itself. 



Figure 2: A friendly path 
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7.1 There exists a friendly path in our cycle 

In this subsection we describe the construction of the canonical path along this 
cycle in the case if a friendly path (going from the main-diagonal to the small- 
diagonal) exists. Let the chords of the existing friendly path correspond to 
the positions ^i, . . . , where Aj = (aj,a|). As we have already mentioned 
we now consider ("comfortably ordered") € x ^ submatrices of our adjacency 
matrices. It is important to remark, that all descriptions of this subsection 
use the language of the matrices Fz (that is the entries can be 0, 1, 2, 3), unless 
it is announced otherwise. 

Before we start we introduce a metric on pairs of positions of this matrix: 
||j4,^|| says how many steps are necessary to go from A to A if in every step 
we can move to a (horizontally or vertically) neighboring position, we cannot 
cross the main-diagonal, finally the position (i, 1) is neighboring to (i,£) and 
analogously (£,«) is neighboring to (1,«). 

By definition our friendly path has the following properties: (i) aj ^ a| 
and aj -f 1 7^ a|. (ii) = 1. Finally (iii) Ai is in distance 1 from the 

main-diagonal, while is in distance 1 from the small-diagonal. 

First we introduce two new structures: 

Definition 7.4. Let 1 < a,/3 < a+1 < (3. We say an £ x ^-matrix M = {rriij) 
is (a, /3)-OK matrix iff 
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(i) ma,fj = 2, 

(ii) TTii.i = for a < i < (3, and ma — 1 for i < a and 13 < 

(iii) rrii^i+i = 1 for a < i < /3, and TOi,i+i = for i < a and (3 <i 

(See Figure Please recall that the entry 2 in Fz is an edge which is missing 
from Z but exists in both G and G' (the off-diagonal entries are the same in 
Mg and Mc)- 



Figure 3: An (q!,/3)-0K matrix and an (/3,q:)-K0 matrix 
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Definition 7.5. Let 1 < a,/3 < £, a+l < (3. We say an £ x ^-matrix M = [rriij) 
is (/3, a)-KO matrix iff 

(i) mp^a = 1, 

(ii) mi^i = for a < i < (3, and ma — 1 for i < a and (3 < i, 

(iii) rrii^i+i = 1 for a < j < /3, and TOi,i+i = for i < a and (3 <i 

(See Figure [3]) Please recall that the entry 1 in Fz is an edge which exists in Z 
but missing from both G and G' . 

Lemma 7.6. Let M = {nii.j) be an (a,l3)-0K matrix and = 3. 

Assume that M' — {m\ j) is an {a — 1, f3 + 2)-OK matrix such that 

(1) = 3 

(2) m-j = m^j ifi^j, i + 1 j, and (ij) ^ (a,^), (a - l,/3 + 2). 

T/ien t/iere exists an absolute constant Q such that one can transform M into 
M' with at most Q swaps. 
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Proof: : It is enough to observe that the symmetric difference of M and M' 
is a cycle which alternates between M and M' . Indeed, in the next figure (1) 
indicates edges in E{M' — M) and (0) indicates edges in E{M — M'). (The non- 
empty positions of this figure are the circled positions in the previous matrix 
M'.) 



E{M' - M) U E{M - M') 






































1 


















1 













1 















































































1 




























































Therefore 

(a - l,a), (a, a), (a,/3,), 
(/?,/?), (/3,/3 + l,),(/? + l,/3-t-l), 
(/3 + l,/3 + 2),(a,/3 + 2) 
is an alternating cycle of length 8 



So the difference of the realizations lay within in a subgraph with 8 vertices 
V . The subgraphs in the two original graphs spanned by V have the same 
(bipartite) degree sequence and they contain alternately the edges of the cycle. 
By Theorem 13.21 we know one can be transformed by swaps into the other one. 
Since the cycle contains four vertices from both classes, and there at most 12 
edges, therefore the number of the necessary swaps (by Corollarv l4.ip is at most 
2x12 therefore 6 = 24 is an upper bound of the number of the necessary swaps. 
□ 

Clearly the same argument gives the following more general lemma. 
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Lemma 7.7. For each natural number u there is a natural number 0„ with 
the following property: assume that M — (mij) is an [a,P)-OK matrix and 
nia'^fji = 3 where 

||(a,/3);(a',/3')|| <u, 
furthermore M' = {'m^ j) is an {a' , f3')-0K matrix such that 

(1) <,0 = 3 

(2) m'^j^mij ifi^j,i + l^j, and {i, j) ^ {a, P), {a',/3'). 
Then at most 0„ swap transforms M into M' . 

Proof: The only difference is that here the symmetric difference of M and M' 
is a cycle of length at most 2 + 2u which alternates between M and M'. □ 

We also have the analogous general result for KO matrices. 

Lemma 7.8. For each natural number u there is a natural number 0^ with 
the following property: assume that M — (rriij) is an (/3,a)-K0 matrix and 
'inpi,a' — where 

II (/3, a); (/?',«') II < u, 
furthermore M' = {m'^ is an {(3' ,a')-KO matrix such that 

(1) = 

(2) rn[j = mij if i ^ j , i + 1 ^ j , and (i,j)^{f3,a), (^',a'). 
Then at most swap transforms M into M' . 

Proof: The proof is very similar to the proof of Lemma [7771 which is left to the 
diligent reader. □ 

Lemma 7.9. Assume that M = i'm-ij) is {a,f3)-0K matrix and mp^2,a-i = 0. 
Assume that M' — {m[ is a (/3 + 2, a — 1)-K0 matrix such that 

(1) <,0 = 3 

(2) m'ij = mij ifi^i, i + 1 j, and ^ {a, I), {(3 + 2, a -I). 

Then there exists a natural number il such that one can transform M into M' 
with at most swaps. 
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Proof: It is enough to observe that the symmetric difference of M and M' is a 
cycle which alternates between M and M' . Indeed, in the next figure values 1 
indicate edges in E{M' — M) and values indicate edges in E{M — M') 



E{M' - M) U E{M - M') 
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Therefore 

(a-l,a-l),(^ + 2,a~l), 

(/3 + 2,/3 + 2),(/3 + l,/3 + 2), 

(/3 + l,/3 + l)(/3,/3 + l), 

(/3,/3),(a,/3), 

(a, a), (a — 1, a) 

is an alternating cycle of 

length 10. 



The proof goes like the proof of Lemma 17.61 The difi'erence of the realizations 
lay within in a subgraph with 10 vertices V . The subgraphs in the two original 
graphs spanned by V have the same (bipartite) degree sequence and they con- 
tain alternately the edges of the cycle. By Theorem 13.21 we know one can be 
transformed by swaps into the other one. Since the cycle contains five vertices 
from both classes, and there are at most 20 edges, the number of the necessary 
swaps (by Corollary ST]) is at most 2 x 20 therefore there exists a constant upper 
bound < 40 of the number of the necessary swaps. □ 

Lemma 7.10. For each natural number u there is a natural number i7„ with 
the following property: assume that M — (rriij) is {a,l3)-0K and mp/^a' = 
where 

||(a,/3);(a',/3')|| < u, 
and M' — {m[ is a {(3' ,a')-KO matrix such that 

(1) <,p = 3 

(2) TOjj = TOij ifi^i, j + 1 ^ j, and (i,j) ^ [a,l), {[3', a'). 
Then at most fl^ swap transforms M into M' . 

Proof: Similar to Lemma [7771 □ 



Now using our friendly path we are going to define a sequence of OK- and KO- 
matrices, such that we can achieve the required edge changes in G getting G' 
along this sequence, using operations described in the previous Lemmas. At 
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first we define a new sequence A[, . . . , A'j^ from Ai, . . . , A^ with tlic following 
way: 

„ _ / A, if Mg(A,) = 0; 

Cousin(A,), if Mg(A,) = 3, 



A'^ 



(7.1) 



where Cousin(A) denote one position where A has a cousin. (So M(3(Cousin(A)) 
Mg(A).) 

Observation 7.11. By definitions, 

(i) ifMciA,) = AfG(A,+i) t/ien A^+J < 3, 

(ii) ifMciA,) ^ MaiA+i) then ||Cousin(AO, J < 3. 



Next we define the matrix sequence Mq 
follows: 



Lq, Li, . . . , L\, La+1 



Mc as 



Definition 7.12. The matrix (i = 1, . . . , A) is defined from the matrix 
by the formulae: 



the (^■)-0K matrix, if L,{A^) = 3; 
the (A'^-KO matrix, if Li{A,) = 0. 



Here all positions {u, v) which are NOT determined by the definitions of the 
OK- and KO-matrices satisfy Li{u,v) = Li-i{u,v). □ 

It is quite clear that (A — 1) consecutive applications of (the appropriate) Lem- 
mas 17.61 - 17.101 will take care the definition of the required swap sub-sequences 
between Li and La. However, the swap-sequence transforming Lq into Li fur- 
thermore the one transforming L/y into Xa+i require special considerations: 

• If Li[Ai) = 3 then there are two possibilities - depending on the position 
of the Cousin(Ai). (The squares denoted with dashed lines contain the 
possible positions of friendly cousins.) 



Case I: 



r — ■ 
1 

1 


3 1 
1 





1 

1 


' 
1 


1 





1 


3 



and Case II. 
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If, however, M{Ai) = then there is only one case: 
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The connecting swap-sequence from the matrix La to La+i (what is Mc) can 
be defined analogously to the previous one. 

Next wo will analyze the behavior of the current matrices M{G, G" , Z) along 
this sub-sequences. At first we consider those Z's which correspond to matrices 
Li. 

Let M be an integer matrix and let M' be a 2 x 2 submatrix of it. If we add 
I's to the values of the positions of one diagonal in M' and — I's to the values of 
the positions of the other diagonal, then the acquired matrix has the same row 
and column sums as M had. Such an operation is called a switch. When our 
matrix M is the adjacency matrix of a degree sequence realization, then any 
swap clearly corresponds to a switch of that matrix. The two matrices are in 
switch-distance 1 from each other. 

Lemma 7.13. All adjacency matrices Li, . . . , L\ are in switch- distance 1 from, 

the adjacency m,atrix of som,e a,ppropriate realizations. 

Proof: We show here the statement for such an Lj where M{Ai) = therefore Lj 
itself is an (A-)-KO-matrix, and where - by definition - Ai = A^ (the other case 
is similar). Due to the dcifinitions Ai originally is not an edge either in G or in 
G' . It belongs to the friendly path, therefore we also know that Cousin(^,) = 
also holds in Mg,Mq. In Li this value is 1, so Ai is an edge in Z. Therefore 
Li which is = Fz looks like the matrix to the left in the following figure (the 
circled element is the cousin of ^j). The corresponding M{G-\-G' — Z) is shown 
on the right hand side: 

M{G + G'-Z) 
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It is clear that adding 1 to the values of the positions A' and Cousin(yl') of 
M{G -\- G' — Z) and subtracting 1 from the other two corners of the spanned 
submatrix constitutes the required switch. □ 
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Lemma 7.14. The realization G can he transformed into the realization G' 
through realizations represented by the adjacency matrices Li (i = \, . . . , K) in 
such a way that the lengths of the swap sub-sequences leading from each Li to 
-Lj+i (where = 1, . . . , Aj can be hounded from above with the absolute constant 
max{03, O3} In the process each arisen matrix M{G + G' — Z) is within a 
constant switch- distance from some vertices in V{G) (that is some realizations 
of the bipartite degree sequence). 

Proof: By Observation l7.1ll anv two A'i and A'i^i are at most distance 3. There- 
fore for each i (where i — 2, . . . , A) the corresponding process chosen among 
Lemma 17.71 Lemma 17.81 and Lemma 17.101 will describe the desired swap sub- 
sequences. The length of any such swap-subsequence is bounded from above by 
max{e3,e(„f}3}. 

Furthermore when in the process the current realization Z corresponds to 
an Li, then Lemma [7 . 1 31 applies . and matrix M{G + G' — Z) has switch-distance 
1 from the adjacency matrix of some realization G V{G). 

Let now Z be an intermediate realization in the process, say, between the 
matrices Li and ii+i: then M{G + G' — Z) can be transform through swaps into 
M(G-\-G' — Lij^i) (assume, this end is the closer one to Z). As we know all swaps 
are specialized switches, and they keep the row and column sums. Combining 
this with the previous paragraph, we have for every Z that M{G + G' — Z) is 
at most i max{83, 63, Vt^} -\- 1 switch distance from some realization e V{G). 
□ 

Key problem 

One can say that we are very close to proving the fast mixing property of our 
Markov process on all bipartite degree sequences: we should prove, that in the 
case when there exists a friendly path from G to G" then for each intermediate 
Z the matrix M{X -\- Y — Z) is in a constant distance from some realization 
G y(G). If we can manage this then we must handle the cases when there are 
no friendly paths. It is somewhat surprising that this second requirement can 
be satisfied successfully (as it will be shown in Subsection 17. 2[) . 

However, we cannot manage to prove the first requirement. The problem is 
the following: we can try to repeat the proof of Lemma !?. 141 but, unfortunately, 
it is not true anymore that for each graph Z, corresponding to a particular 
matrix Li, the matrix M{X -\-Y — Z) is also in distance 1 from some realization 
e V{G). 

In the realizations G and G" all chord have the same values, but this is not 
the case for realizations X and Y. The edges in E{X — Y)\J E{Y — X) belong 
to only one of them. Therefore if a swap turns an entry to 2 in M{G + G' — Z) 
then this entry originally was 1: the edge belonged to G and G' and Z as well. 
Therefore its cousin bears the entry 1 (also belonged to G and G' and Z as 
well). So this entry was appropriate to perform a switch to turn the matrix 
under investigation into the adjacency matrix of a realization. 
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However, if the cousin entry is in M{X + Y — Z) (this edge belongs only 
to one of realizations X and Y, say, it belongs to X only), then the required 
switch cannot be performed. (The value —1 can be cause a similar problem 
and can be handled similarly as this case.) 

A good solution for this particular problem would probably ends up in a 
complete proof of the fast mixing property. 



The following observation is enough to handle the switch-distance problem for 
M{X + y — Z) in semi-regular bipartite degree sequences: recall, that (a, b) is 
semi-regular if in a all degrees are the same, while in b can be anything. 

Lemma 7.15. Assume that our bipartite degree sequence (a, b) is semi-regular. 
Then the statement of Lemma Vl applies for the matrices M[X + Y — Z) as 
well. 

Proof: We follow the proof of Lemma 17.141 To do so the only requirement is 
to show (somewhat loosely) that the matrices M{X -\-Y — Li) are in a constant 
switch-distance from the adjacency matrix of some realizations. As we know 
any of these matrices contains exactly one entry of value different from 1 and 
0. So consider a particular Li and assume that this "extra" value in this case 
is a 2. If the switch, described in the proof of Lemma [7.131 is also a possible 
switch in AI{X -\-Y — Z) then we are ready. If this not the case then the entry 
(with value 1 in matrix M{G -\- G' - Li)) has value in M{X + Y - L,). (In 
this case, as we discussed it previously, the corresponding edge is missing from 
Y.) Let this corresponding edge be {u, v), then this entry in M{X -\-Y — Li) is 

0. Since the column sums are fixed in these matrices, they are the same (and 
equal to entries in a). 

Now vertex v has degree at least 2 (it is a vertex on cycle C and it also end 
point of at least one chord of C in X). Therefore the raw v contains some Is. 
One of them is {w,v) (this w cannot be the column of the 2, since the entry 
there is due that it belongs to the originally intended switch). Now by the 
pigeonhole principle there is a raw z such that M{w,z) = and M(u, z) — 1. 
Therefore the u, w; v, z switch (actually this is a swap) will change M(u, v) into 

1, and now the original switch finishes the job. The matrix M{X -\-Y — Li) is 
in switch-distance at most 2 from the adjacency matrix of some realization. □ 

7.2 Their is no friendly path in our cycle 

In this subsection we discuss the case where there is no friendly path in Mq. 
Our plan is this: at first we show that the non-existence of the friendly paths 
yields a strong structural property of the matrix Mq- Using this property we 
can divide our problem into two smaller ones, where one of the smaller matrices 
possesses a suitable friendly path. So we can solve our original problem in a 
recursive manner. 

This recursive approach must be carried out with caution: a careless " greedy" 
algorithm can increase the switch-distances very fast. We will deal with this 
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problem using a simple "fine tuning" (which is described at the end of this 
subsection). 

We start with some further notions and notations. 

Definition 7.16. Let M be an ^ x ^ matrix, then the positions + 

1, i — 1), (i + 2, i — 2), . . . , + 1) form the ith down-line of the matrix. 
(The arithmetic operations are thought to be considered modulo ^, that is, for 
example, £ + I ~ I while 1 — 3 = ^ — 2. The positions {i, i), {i — l,i + 1), {i — 

2, j + 2), . . . , (j', j' — 1) form the ith up-line of the matrix. 

The lines have orientations: the down-lines connect the main-diagonal to 
the small-diagonal, while an up-line connects the small-diagonal to the main- 
diagonal. 

Definition 7.17. A set T of positions of an £ x £ matrix called connected if 
a chess king, staying inside T, can visit all elements of T. 

The following lemma is a well-known version of the classical Steinhaus lemma 
(see Hg). 

Lemma 7.18. Assume that the positions of an £ x £ matrix M are colored for 
white and black. Then either the king has a white path from the main diagonal 
to the small diagonal, or there is a connected set T of black positions which 
intersects all rook's path from the main diagonal to the small diagonal. □ 

We use the previous result without proof. The set T, which was identified 
in the previous lemma, will be called a Steinhaus set. 

Definition 7.19. The cousin-set £(u, v) is the set of the off diagonal cousins 
of the position {u,v). If T is a set of positions, then the cousin set £(T) is 
defined as U{e(e) : e G T}. 

Lemma 7.20. Assume that in the matrix Mq there is a connected set T of 
non-friendly positions. Then the type of all positions in the cousin-set €{T) are 
the same {all chords have the same type). All edges in T have the opposite type. 

Proof: W.l.o.g. we may assume that a position P in T has type 0, then all 
positions in its cousin-set must have type 1. However, for each other position 
P' in r, which can be reached from P in one king step, the cousin sets £(P) 
and C{P') have common position(s). Therefore all type in those two cousin sets 
must be the same (1), therefore both positions P and P' have the same type 
(0) as weU. □ 

Lemma 7.21. Let T be a Steinhaus set of unfriendly positions in Mc, then its 
cousin-set £(T) intersects all down-lines and up-lines. 

Proof: Actually we can prove more: namely that any king-path from the main- 
diagonal to the small-diagonal intersects the cousin-set £(r). We will argue by 
contradiction, assume that there exists a king-path which does not cross the 
cousin-set. For position P let denote I{P) be those positions whose cousin-sets 
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contain P. Then if P e king path, then I{P) is disjoint from C(T) (otherwise 
the king-path intersects T). But any two neighboring positions in the king- 
path have overlapping I{P) sets. Therefore we have an overlapping sequence 
of I{P) sets, connecting the main-diagonal to the small-diagonal. And such an 
overlapping sequence clearly contains a rock-path, connecting the main-diagonal 
to the small-diagonal. A contradiction. Finally it is clear, that every down- and 
up-line forms a required king-path. □ 

Lemma 7.22. We assume that in the matrix Mq there is no friendly path. 
Then for each i (i ~ 1, . . . ,£) there exists a pair of s.t. < j' — j < 1 and 
all entries {i + l,i — 1, . . . , {i + j,i — j) have the same type t, furthermore the 
entries (i — 1, i + 1) , . . . , (i — j' -\- I) , i + j' ^ 1) have type 1 — t while {i + j' ,i — j') 
has type t again. 

Proof: Assume indirectly, there is no such j for a particular i. W.O.L.G. we 
may assume, that M{i -I- 1, j — 1) = 0. Then, by the assumption, M{i — 1, i -I- 1) = 
M{i — 2,i + 2) = 1 must hold. Then, again by our assumption, M{i + 2,i — 2) = 
must hold, etc. All entries along the down-line are 0, while all entries along the 
up-line must be 1. However both lines intersect (see Lemma l7.2ip the cousin-set 
£(T) of the Steinhaus set T. But, by Lemma [7.201 all its entries have the same 
type. A contradiction. □ 

Corollary 7.23. If conditions of Lemma \ 7. 2'2\ hold, and j' > 2 for a particular 
i (in which case j' — j the submatrix spanned by {i+j, i—j) and {i—j, i+j) 
contains friendly path(s). 

Proof: We argue by contradiction: assume that the submatrix does not contain 
a friendly path. Then - due to Lemma 17.181 - it contains a Steinhaus set. Due 
to Lemma [7.201 in its cousin-set £(T) - which intersects all down- and up-lines 
- all positions have the same type. But it contradicts to the fact, that in the ith 
down- line all positions have type t, while in the ith up-line all positions have 
type 1 — t. A contradiction, again. □ 

That finishes the preliminaries what are needed to describe our recursive algo- 
rithm, which is essentially a divide and conquer approach. Due to the previous 
fact here we should handle separately two possibilities: when j' — 1 and when 
/ > 2. We start with the 

First possibility: assume that for a particular i our j' = 1. We should take 
care for two cases: 

Case 1: If Mq = (i -I- l,j — 1) = Mcii — l,i -I- 1) — 1 then we are in an 
easily handleable situation: at first we swap the quartet Mi, Mi+i; Vi, w^+i. (The 
dash-ed square in our illustration. Here we use the adjacency matrix Mg-) The 
entries [i — l,i),(i,i) and (i,j -I- 1) have the required types. However, entry 
M(« — 1,1 + 1) ~0 therefore along the procedure we should take care to change 
it back to its original value. 
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The remaining subproblem, indicated by thick black hnes (for a formal de- 
scription see case j' > 2), fortmiately is already in the required form. Indeed: its 
main-diagonal contains only Is, while its small-diagonal is full with Os. Denote 
the alternating cycle of this smaller problem by C". The (recursive) solution of 
the subproblem C will switch the value of M(i — 1, i -I- 1) automatically back 
to 1. Since the recursive procedure can use any down- and up-lines (see Corol- 
lary [7jl3]), therefore we can take care, that this switch-back will happen in the 
next recursion. 

It is important to recognize that matrices AI{G + G' — Z) and M{X + Y — Z) 
may have contain 2 at the position (i — 1, i -|- 1). Fortunately this "problematic" 
entry will be present only along one recursive step. Furthermore this entry will 
increase the switch-distance of the current M with at most one: the positions 
(i — and + 1) (outside of our subproblem), provides a suitable 

switch to handle the entry 2 at position (i — 1, i -I- 1) . 

Case 2: Now we have Mq = {i + l,i — 1) = Mc{i — 1, i -I- 1) =0. Here we 
perform two swaps, as it is shown below (the places of the swaps are denoted 
with dash-ed squares): 
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If Ad{X + Y — Z^{i + 1, i — 1) = — 1 holds, then it increases the switch-distance 
of the current M with at most one (since it can be directly back-swapped). 
The result of the second swap (after which the previous problem is just solved 
automatically), together with our further strategy is shown below: 
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Here we distinguish between two 
cases, according to the value 

Mcii -2,i + 2)=x. 

This value can he x — 1 and a; = 0. 



In the case of a; = 1 we perform a second swap, which results in a subproblem 
with a friendly path (the swap shown on the left side of Figure |4j while the right 
hand side indicates the two new subproblems): 



Figure 4: The case of a: = 1 




















1 

















1 

















1 








] 1 






1 


1 















1 


1 














1 

















1° 







1 1 













1 
















1 




































1 

















1 

















1 















1 


1 












1 


1 


] 














1 

1 1 
1 

















1 


1 
1- 




1 













1 
















1 



















On the RHS of Figure SI the dash-ed submatrix is in the right form: the main 
diagonal contains all I's, while each entry in the small-diagonal is 0. The second 
subproblem (indicated with the thick black lines, the four pieces fit together to 
a square matrix, again in the right form. We will process the second ("dash- 
ed") subproblem along the up-line, containing position {i — 2,i + 2), so the 
only currently improper entry will have the right value at the end of the next 
recursion step (that is it will be swapped back to its original value). (Here we 
used again CoroUarv l7.23n 
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As it happened before M{X + Y — Z) may contain 2 at the position {i — 
2,i + 2). Again this increases the switch-distance with at most one, since the 
positions {i — 2,i— 1), + and {i + l,i + 2) are not in our subproblem. 



Finally it can happen, that x = 0. Then we can define the following subproblem: 
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This figure shows the new sub- 
problem (indicating with thick 
black lines) is in the right form 
again. We will process the sub- 
problem along the up-line, con- 
taining position (i — 2, i + 2) (so 
the only currently improper entry 
will have the right value at the end 
of the next recursion step). 



Here, again, we may confront the fact, that M(X + 
Y - Z){i + l,i - 1) = -1. Then we should con- 
sider the alternating cycle shown in the figure. All 
elements of the cycle, except {i + l,i — 1), is in 
the main- and small-diagonal, therefore along this 
cycle wc can swap that entry into range within a 
small number (say 5) of steps. This will increase the 
switch-distance of M with at most 5. 

Wc nm the first recursion on the subproblem along the ith up-linc, therefore 
the sub-subproblem with friendly path will contain the position (i — j' ,i + j'). 
Therefore when we finish the first recursion, our adjacency matrix Mq will be 
in the following form: (the figure on the left): 



3 














1 







2 
















1 







1 









1 
















1 


2 


3 

















1 

















1 

















1 











1 






















1 


1 
















1 

1 1 
1 





1 
1 















1 
1 




r 

1 1 










1 





















































1 

















1 

















1 











] 1 






1 


















1 


1 














1 

















1° 







1 1 










1 







































We have seen how one can handle the switch-distance of our matrix, if position 
{i + l,i — 1) is problematic (in M it is a —1 there) but position {i — 2,i + 2) 
is O.K. On the other, if M{X + Y - Z)){i -2,i + 2) = -1 then the swap on 
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the positions (« — 2, j — 1), (i + 1, « + 2); (i + 1, i — 1), (i — 2, i + 2) change both 
(i + 1, i — 1) and (i — 2, i + 2^ into 0. For {i + l,i — 1) that was the original type 
- so it cannot be wrong in M. 



After that we perform the swap on the positions {i — 

2, 1 - 1), (? + 1, 1 + 2); (i + 1, i - 1), (i - 2, 1 + 2) (these 
are the corners of the dash-ed square in the figure on 
the upper right). The result is shown to the right: 
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This completes our handling on the First Possibility, that is when for our i we 
have the value j' = 1. Now we turn to the other (and probably more common) 

configiiration: 

Second Possibility: Wc have j' > 2. Unfortunately, the situation can be 
more complicated in this case duo to the switch-distance of M. We overcome 
this problem by showing at first the general structiire of the process, and later 
we give the necessary fine-tuning to ensure the low switch-distance. 

In our current alternating cycle (lying in the symmetric difference of G and G") 
there is no friendly path, therefore there is a T Steinhaus set in the subgraph of 
Z, spanned by the vertices of the cycle. Now fix a particular i and assume that 
the j' belongs to this z is > 2. We should distinguish between two cases: where 
the down-line start with the value t = 1 or with t = 0. 

Case 1: t = 1 The first figure below shows the structure of matrix Mg- The 

dash-cd square is the first subproblem to deal with, while the thick black lines 
indicate the second subproblem. However, before we start the processing the 
subproblems, we have to perform a swap. The corners of the thin black square 
shows the positions of the swap. After that, the first subproblem (indicated by 
the dash-ed square) is in the right form. (See the figure on the right.) 



















n 


1 
























1 































1 





















1 


















' 


1 




1 












( 


)- 


















j-^ 















1 

















1 
















1 


































1 


1 








1 ° 









1 ' 


















1 

















1 














I 


1 




1 










1 


[ 1 








! 









1 

















1 
















1 



















1 1 
Finishing the first subproblem, we have the following Mc adjacency matrix: 
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As it can be seen, after the first 
phase, all entries in the dash-ed 
are in their required types: the 
small-diagonal consists of Is (in- 
cluding position (i + 2, i — 2) which 
in that way is back to its origi- 
nal type), while the main-diagonal 
consists of only Os. 
The second subproblem (indicated 
by the thick black lines) in the 
right form now (including position 
(i — 3; ?' + 3) which is sitting on the 
^ small-diagonal). 
After completing the solution of the black subproblem, all entries in the matrix 
will be in exactly the required type. We start processing the black subproblem 
on the ith up-line, therefore the actual types of positions (« — 3, i -h 3) and 
{i + 2,i — 2) can be described as follows: Position {i + 2,i — 2) has opposite type 
after the very first swap, then while processing the dash-cd subproblem it may 
changes between and 1. Finishing the dash-ed subproblem, it will be in the 
type as it starts. 

Position (i — 3, i + 3) will be in type all the way in the dash-cd phase, while 
within the black phase it will change between 1 and 0. At the end, as we already 
mentioned, it 1. 
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Case 2: t = The first figure below shows the structure of matrix Mq- The 
dash-ed square is the first subproblem to deal with, while the thick black lines 

indicate the second subproblem. They can process without any preprocessing. 
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At the end everything will be in the right type, except the four position, showed 
by the thin black square (below, left side). We can finish the process with that 
swap. 
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While the overall structure of our plan is clear, we may meet problems along 
this procedure. Its reason is that we must be able to control the switch-distance 
of our M{X + Y — Z) (we will use here simply M) from the adjacency matrix of 
some realization. There are two neuralgic points: both the positions {i+j, i—j) 
and {i — + j') niay contain —1, or both may contain 2. When we start a 
new subproblem, then their types always provide suitable switch for the control 
(as it was seen before). However, when we proceed along our subproblem, then 
it can happen that one of the problematic position changes its value, while the 
other does not. But in this case the switching which was previously available are 
not useable anymore. Next we describe how we can fine tuning our procedure 
to avoid this trap. 

As we know the first subproblem contains a friendly path, and for easier 
reference let call its problematic position Pi. We also know that second sub- 
problem contains a problematic position, P2, and probably we have to divide 
this subproblem into two smaller ones. If so, then the first of them becomes the 
new second subproblem, which contains P2 and possesses a friendly path, while 
the third subproblem contains another problematic position, P3. 

Fine tuning: 

1. We begin our swap sequence along the first subproblem but we stop just 
before we face the swap which changes the value of Pi . 

2. Next we continue with the swap sequence of the second problem and we 
stop before we should perform a swap on P2. 

3. Now we finish the swap sequence of the first subproblem. 

4. After that we focus on the second subproblem. Dealing effectively with 
this, we need to prepare the third subproblem similarly as we did with the 
second one, when we were working on the first one. Therefore we begin 
the swap sequence of the third subproblem but we stop it before the first 
swap would be carried out on P3. 
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5. And if now we just rename our two active subproblems as first and sec- 
ond subproblem, we are back to a situation, which is equivalent to the 
beginning of the third stage. 

Doing this refined algorithm we have the opportunity that at every particular 
moment we can control with a small number of swaps the M values at positions 



The size of the parameter set B: In this paper we proposed the friendly 
path method to construct suitable canonical path systems for the application 
of the Sinclair's method. Lemma 16.31 postulated several assumptions which 
were managed throughout the construction. However we did not check the 
assumption about the size of the parameter set B yet. Here we supply this last 
missing part. 

Lemma 7.24. The size of the parameter set B is poly(n). 
Proof: The parameter set contains the following parameters: 

1. The integer j describing the ordinal number of the current step. Since no 
path is longer than Cn^ therefore 1 < j < also holds. 

2. The edge pair a, b describing the segment [a, b] of the current alternating 
cycle which is already switched from E{G) toE{G'). Since n is a trivial 
upper bound of the length of any cycle, therefore there are at most O(n^) 
possible edge pairs. 

3. We should list those chords of the current alternate cycle which have 
different types in X and Y and touched by the construction. At the 
very first and last steps of the construction of the swap sequence along 
a friendly path there may be at most one "bad" chord. While along the 
steps described in Lemmas |7. 6117. 7|[7?51 there are at most max{6,6„,6^} 
bad chords. (In each application of those sub-procedures at most 2 bad 
chords can come into existence and the same amount can be omitted. At 
the end of those applications there are always at most one bad chord. 

4. Finally at the fine tuning a constant number C of bad chords may come 



All together this gives an poly(n) upper bound on the size of B. (Here we gave 
only very rough upper bounds but they were sufficient for our purposes.) □ 

8 Conclusions 

In this paper we introduced a new way to generate canonical path systems for 
the Sinclair's method to be applied for the Markov chain determined by a given 
bipartite degree sequence (defined by Kannan, Tetali and Vempala [5]). In the 
process no specific assumptions on the degree sequence are adopted, but showing 
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the validity of the argument requires the assumption of semi- regularity (at least 
for now) . See Lemma 17.151 and the Key Problem. However we think that this 
technical limitation can be eased by the mean some new, but not necessarily 
breath-taking ideas. Therefore here we summarize the most important feature 
of our " friendly path system" method: 

As Lemma 17.141 and the procedure described in Subsection 17.21 proves the 
following statement holds: 

Theorem 8.1. Let G and G' be two realizations of an arbitrary bipartite 
degree sequence, such that the edges E{GAG') form one alternating cycle. Then 
realization G can be transformed into the realization G' through realizations in 
such a way that the length of the swap sequence leading from G to G' can be 
bounded from above with some absolute constant. In the process each arisen 
matrix M{G + G' — Z) is within a constant switch- distance from some vertices 
in V{G) (that is some realizations of the bipartite degree sequence). 
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